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ABSTRACT 

Utilizing the CLASS statistical sample, we investigate the constraint of the 
splitting angle statistic of strong gravitational lenses(SGL) on the equation-of- 
state parameter w = p/p of the dark energy in the flat cold dark matter cosmol- 
ogy. Through the comoving number density of dark halos described by Press- 
Schechter theory, dark energy affects the efficiency with which dark-matter con- 
centrations produce strong lensing signals. The constraints on both constant 
w and time-varying w(z) = w + w a z/(l + z) from the SGL splitting angle 
statistic are consistently obtained by adopting a two model combined mecha- 
nism of dark halo density profile matched at the mass scale M c . Our main 
observations are: (a) the resulting model parameter M c is found to be M c ~ 1.4 
for both constant w and time- varying w(z), which is larger than M c ~ 1 ob- 
tained in literatures; (b) the fitting results for the constant w are found to be 
w = — 0. 891^26 an d w — ~ 0-94lo!i6 for the source redshift distributions of the 
Gaussian models g(z s ) and g c (z s ) respectively, which are consistent with the 
ACDM at 95% C.L; (c) the time-varying w(z) is found to be for cr 8 = 0.74: 
{M c ;w ,w a ) = (1.36; -0.92, -1.31) and (M c ; w Q , w a ) = (1.38; -0.89, -1.21) for 
g(z s ) and g c (z s ) respectively, the influence of a s is investigated and found to 
be sizable for o"g = 0.74 ~ 0.90. After marginalizing the likelihood functions 
over the cosmological parameters (fl,M,h,as) and the model parameter M c , we 
find that the data of SGL splitting angle statistic lead to the best fit results 
{w ,w a ) = (-0.88^,-1.551^) and (w ,w a ) = (-0.91±?;Sg, -l-60±|;i?) for 
g(z s ) and g c (z s ) respectively. 



Subject headings: cosmological parameters — cosmology: observations — cosmology: 
theory — gravitational lensing — quasars: general 
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1. INTRODUCTION 

Since the direct confirmation of the presence of dark energy by Type la supernovae 
(SNe la) observation (Riess et.al 1998), the investigation of its property is one of the most 
important object in cosmology. Many theoretical models have been developed to explain or 
describe dark energy, which is widely believed as the main component of the cosmological 
energy today. By now one key to the question seems to be the precise measurement to 
the equation-of-state parameter w = p/p of the dark energy(see a review paper, Peebles & 
Ratra 2002). Compared with the data, the limit to the value of parameter w is continuously 
improved by many experimental groups, including the type la supernova(Riess et al. 2004), 
the Cosmic Microwave Background(CMB)(spergel et al. 2006), and the weak gravitational 
lenses(WGL) (Weinberg & Kamionkowski 2002). As a complement to these observations, we 
shall utilize the splitting angle statistic of strong gravitational lenses(SGL) to quantification- 
ally investigate its constraint on the parameter w. 

Light lines traversing in the universe are attracted and refracted by the gravitational 
force of the galaxies on its path, which bring us the signal of the SGL effect, one of which 
is the multiple images of a single far galaxy. Through comparing the observed number of 
lenses with the theoretical expected result as a function of image separation and cosmological 
parameters, it enables us to determine the allowed range of the parameter w. Linder (2004) 
demonstrated that with the addition of strong lensing image separation measurements, the 
estimates for time-varying w(z) = wo + w a z/(l + z) from CMB, SNe la, and WGL could 
be improved modestly. To see that, we shall carefully in this note investigate the analytic 
process and the power of SGL data only. There are numerous works studying on the relation 
between DE and SGL splitting angle statistic (e.g., Porciani & Madau 2000; Li & Ostriker 
2002; Kuhlen, Keeton & Madau 2004), but the exact mass density profile of dark halos is 
unknown yet, which produces the most theoretical uncertainty of the analytic process. The 
two most widely-used density profiles are the singular isothermal sphere (SIS) profile and 
the Navarro-Frenk- White (NFW) profile. Some analyses were based solely on the SIS model 
(e.g., Chae et al. 2002). While from the analyzes by Li & Ostriker (2002, 2003), and also by 
Sarbu, Rusin & Ma (2001), it gave a convictive illumination that a combined mechanism of 
at least two models can effectively reproduce the observed curve of the lensing probability 
P{> 0) to the image splitting angle 9. To achieve this, an additional parameter M c is 
introduced to divide the mass scale of dark halos into different parts as a certain density 
profile is thought to work only on each part. 

The rate of structure growth, which determines the number density of dark halos as the 
SGL lenses, is very sensitive to the normalization parameter of the matter power spectrum, 
a 8 . The value of cr 8 is often related to the matter density Qm and constrained by the CMB or 
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cluster abundance observation. Several years ago, in the ACDM universe with Q M = 0.27, 
a 8 was found to be larger than 0.9 (e.g., Spergel et al. 2003; Wang & Steinhardt 1998). 
Then the data of Wilkinson Microwave Anisotropy Probe (WMAP) Three Year (Spergel 
et al. 2007) provided a noticeably smaller value a 8 = 0.74 ± 0.05 in comparison with the 
first-year data erg = 0.92 ± 0.10. By utilizing the recent clustering results of XMM- Newton 
soft (0.5-2 keV) X-ray sources, the X-ray clustering, and SNIa data, Basilakos and Plionis 
(2007) showed that the a s 0.73 for Vl M = 0.26 and w = —0.90, which is consistent with 
the new result. A smaller cr 8 implies less structure growth at late times and less lensing 
probability for a single source. 

In this paper, by using the CLASS statistical sample (Browne et al.2003), we are go- 
ing to show how the constraint on the dark energy equation of state parameters w(z) can 
reliably be obtained for the cases with the assumption of constant w and the time-varying 
parameterization w(z) = wo + w a z/(l + z), respectively. We highlight three issues which 
have not previously been investigated. First, we investigate the influences of parameter w 
on the every step of the lensing process to find the most important point where the change 
of parameter w shows its effect. Second, by comparing the results for different distributions 
of sources redshift which have been used in previous works(e.g., Chae et al. 2002; Li & 
Ostriker 2002), we illustrate the quantitative influence that is introduced by the uncertainty 
of the source distribution. Third, we focus on the constraint of the data on the time-varying 
parameterization w(z). Our paper is organized as follows: Sect. 2 outlines the cosmo logical 
model, the mass fluctuations, and the Press- Schechter function used in our calculation. Sect. 
3 describes the density profiles and the lensing probabilities. Sect. 4 gives our data analysis 
and numerical results. The conclusions are presented in the last section. 

2. BASIC CONSIDERATIONS 

In this section we shall describe some ingredients used in our calculations. 

2.1. Cosmological Model and Mass Fluctuations 

During the last two decades more and more observational evidence suggests that our 
universe at present is accelerated expanding and dominated by a spatially smooth component 
with negative pressure, so called dark energy. Besides the common cosmological constant 
supposition, an attractive alternative candidate for dark energy is the potential energy of 
a slow-varying scalar field, which is conveniently parameterized through w — p/p. The 
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conventional scalar field models, i.e. quintessence models, with a positive kinetic energy 
term in the field Lagrangian require w > — 1 (e.g., Ratra & Peebles 1988; Caldwell, Dave, 
& Steinhardt 1998), and phantom dark energy , adopting alternatively a negative kinetic 
energy term, gives the parameter space w < — 1 (e.g., Caldwell 2002; Carroll, Hoffmann, & 
Trodden 2003; Cline, Jeon, & Moore 2003 ). To fit the data of the observation and give 
the allowed parameter space, we use two typical parameterizations, i.e., constant w and 
time- varying w(z), as follows: 



The simple form w(z) = wo + wi z is not favored because the value of w(z) runs to infinite 
when redshift z goes to be infinite. Throughout this paper, we assume a flat cold dark matter 
universe with present matter density relative to critical density VLm = 0.24 and the Hubble 
parameter h = 0.73 (Spergel et al. 2007). 

Dark halos of galaxies and galaxy clusters are formed through linear growth and non- 
linear collapse of primordial fluctuations of matter in the early universe. The time-varying 
state-of-equation parameter w(z) is usually realized by models with slow-varying scalar fields, 
while these scalar fields of dark energy begin to cluster gravitationally and contribute to the 
perturbation spectrum only at the very large spatial scale L > lOOMpc, which corresponds 
to a very small wavenumber k < O.OlMpc -1 (e.g., Ma et al. 1999). In our present consid- 
erations, the most concerned length scale I < IMpc is far less than the length scale L and 
the concerned fluctuations are obtained by integrating the mass power spectrum k 2 P{k) up 
to a larger value of the wavenumber k > l.OMpc -1 . Noticing the fact that the increasing 
function k 2 P(k) oc k 3 when k « l.OMpc -1 , so the contributions from the integral region 
k < O.OlMpc -1 are very small, namely the contribution of the dark energy cluster to mass 
fluctuations can be neglected. Thus we can directly utilize the mass power spectrum of 
ACDM cosmology. In our calculation, the linear CDM power spectrum is computed by 
adopting the fitting formulae given by (Eisenstein & Hu 1999) 



The initial power spectrum index n is fixed to be n = 1. T is the transfer function 



Casel : w = constant 



Casell : w(z) = w + w a z/(l + z), w, 
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dw(z)/dz\ z=0 . 
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The parameter 5h is the amplitude of perturbations on the horizon scale today and 
related to the rms density fluctuations in spheres of radius r 8 = 8/i _1 Mpc, so called a 8 by: 
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. To show the 



where W(kr) is the top-hat window function: W(kr) = 3 
power of SGL data only, we don't use any analytic fitted form of the parameter cig, expressed 
by Qri and w. Alternatively , unless special clarification we choose to normalize the power 
spectrum to a§ = 0.74, the best-fit value given by WMAP Three Year data (spergel et al. 
2006). 

The linear growth function T>(z) is proportional to the linear density perturbation 5 = 
Spm/pm- The evolution of linear perturbation is: 



6 + 2-6 = 4:7iGp M 6 



(6) 



where a is the scale factor a = (1 + z)~ , dot means derivative with respect to physical time 
t, the background matter density p M = p (l + z) 3 , p = ^A/Pcrit,o an d p cr it,o = 3Hq/(8ttG) is 
today's critical mass density in the universe. Then with the definition D(z) = 6(z)/6(z = 0), 
we can obtain the equation of T>(a): 



d 2 V 3 tt^ 3 dV. . W1 

— — = — -v — i - w(a)(i - n)\ 

da 2 2 a 2 2a da 1 V A n 
where Q is the matter density parameter Q = Qm(1 + z) 3 /(H/Hq) 2 and 



(7) 



H 



aHn 



'n M (l + z) 3 + n DE Exp (J dz'3[l + w(z'))/(l + z') 



Here Hq is the present Hubble constant Hq = 100/ikms^ 1 Mpc _1 and Qde = (1 — for a 
flat universe. With two boundary conditions: V(a)\ a=0 = and V(a)\ a=1 = 1, Equation (|7j) 
can be calculated out numerically. The second boundary condition means T>(a) is normalized 
to Via = !) = !. 
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2.2. Mass Function and Spherical Collapse Approximation 



According to the Press-Schechter theory, the comoving number density of dark halos 
virialized by redshift z with mass in the range (M, M + dM) is given by 



n(M,z)dM = ^f(M,z)dM. 



(8) 



f(M, z) is the Press-Schechter function. We utilize the modified form by Sheth & Tormen 
(1999) 



f(M,z) = 



0.383 5 r dA 



A 2 dM 



A 5 



0.7075 2 



0.3 



x exp 



0.707 



f°° dk 

A 2 (M,z) = J T A k (k,z)W 2 (kr) 



(9) 



(10) 



where A is the variance of the fluctuations in a sphere containing a mean mass M, and M 
is related to the length scale r via 

M = yr 3 pM- 

The parameter 5 c (z) is the linear overdensity threshold for a spherical collapse by redshift 
z. The matter with the overdensity in a certain scale of the early universe would undergo 
the density growth and the spatial scale reducing. When its average matter density reaches 
S c (z), virialization starts and then a dark halo is formed. In this paper we follow Wang & 
Steinhardt (1998) and Weinberg & Kamionkowski (2003) to calculate the S c (z). Under the 
approximation of spherical tophat collapse and labeling R as the spatial length scale for a 
halo with a certain mass, the collapse process is determined by the Friedmann equation 

2 



8ttG 



(Pm + Pde) , 



and the time-time component of the Einstein equations 
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(11) 



(12) 



where Pde is the energy density of dark energy and phaio is the uniform matter density in the 
scale r . When the scale factor a(z) is very small, i.e., a(z) — > or z ~ zq — > oo, the equivalent 
linear overdensity threshold 5 c (z) at z can approximately be evaluated through the dark 
halo density Phaio an d the background matter density pm, i-e., S c (z ) ~ (phaio/PM — 1)- 
Then utilizing the two boundary conditions: r(a)\ a =o = and dr/da\ a=ata = with a ta the 
scale factor at the turn- around time, the function S c (z) can be calculated out numerically as 
follows 

$c(z) = &c{zq)D(z)/D(z ), 5 C (Z ) ~ (phalo/pM ~ 1) 
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Note that the turn-around time tta is determined through the virial time t V { r when the 
overdensity matter starts to form dark halos, i.e., t ta = t vir /2, which is corresponding to virial 
redshift z according to the Press-Schechter theory. In this sense, the boundary condition 
dr I 'da\ a=ata = is related to the virial redshift z, so the resulting scale r(a) of dark halos 
will depend on the virial redshift z. As a consequence, the dark halo density phaio actually 
relies on the virial redshift z. 

Denote A vir = Phaiol Pm the ratio of the cluster to the background density. Then the 
nonlinear overdensity A vir (z) can also be calculated out directly from the Eq. llTTT) and Eq. (IT21) 
with their two boundary conditions and the 5 c (zo) given above. 

In Fig. 1, we plot the Press-Schechter function / against the mass M(10 15 h~ l M Q ) 
of dark halos at redshift z = 0.0, 1.5, and 3.0, respectively. The dash, solid, dot, and 
dash dot curves are for the typical cases of w — —0.5, w = —1.0, w = —1.5, and w(z) = 
—1.0 — z/{l + z). Fig. 2 shows the P-S function / as a function of the redshift z with the 
dark halo mass M(10 15 h~ 1 M Q ) = 0.001, 0.1, and 10 , respectively. Notice the differences of 
the ordinate length scales. From the two figures, we can see that for larger M and smaller 
z, the P-S function / is more sensitive to the change of the parameter w. The explanation is 
clear: (1) more time passed for the fluctuations to evolve and form dark halos of larger mass, 
so larger M corresponds to later time; (2) in the early time, the relative proportions of dark 
energy in the total energy of universe was small and so is its influence on the spatial-time 
geometry. We can see that even for a small M ~ 0.001, the change of w can bring somewhat 
significant shift of function / in the most range of redshift z. This produces the most power 
of SGL data to constrain the parameter w. 



3. DENSITY PROFILE AND LENSING PROBABILITY 

The SGL lensing efficiency is very sensitive to the density profile of dark halos: under 
the same conditions, the efficiency of SIS model is larger, at least by one order of magnitude, 
than that of NFW model for the image separation angle 58 < 30". In this section, we shall 
discuss the influence of different w on the lensing probability P{> A#o) for SIS profile and 
NFW profile, respectively. 
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3.1. SIS Profile as a Lens 

The SIS profile has a simple spherically symmetric form (Schneider, Ehlers, & Falco 



M=£bh- (13) 



where a v is the velocity dispersion, which can be related to the mass M of a dark halo via 
al = GM/2 r vir after integrating the density function p(r) from r = to r = r vir . Here r vir 
is the viral radius of a dark halo, which is commonly defined by demanding that the mean 
density within the virial radius of the halo be a factor A vir times larger than the background 
density, pm, i.e. M = ^fA vir p cr r1 ir . Eliminating the dependance on the r vir , we get 

al( 6 N 1/2 



M=-£\ t " . (14) 

G \7rA vir G p cr J 

This profile is supported by the observed flat rotation curves of the spiral galaxies and 
is widely utilized in the gravitational lensing anslysis. Due to its symmetry, the lensing 
analysis is quite easy. Integrate the density component along the line of sight and then we 
get its surface mass density 

m = £\. ("> 

where £ = |£| and £ is the position vector in the lens plane. The lensing equation is given 
by 

ri = - Di s a{Z) , (16) 

where r] is the source position. Dg, , and D^ s are the angular-diameter distances from 
the observer to the source, from the observer to the lens and from the lens to the source. 
The angle <*(£) is the gravitational deflection angle. For a circularly-symmetric surface 
mass density, £(£)= £(£), images appear on the plane defined by the observer point, the 
lens center, and the source position, and the angle a(£) is given by 

«(o = ^? fmodt (17) 

c s Jo 

To simplify the lensing equation, we define the length scales in the lens plane and the source 
plane as 



(7, 



-sf- * = «»4' (18) 
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Then the position vector of a point in the lens plane or source plane is £ = x£o or 77 = yr/o. 
After the reduction, the lensing equation for an SIS lens is given by 

y = X -W. (19) 

X 

It is easy to see that when \y\ < 1, i.e. \x\ < 1, a single source has double images with the 
separation Ax = 2 and the splitting angle 

AO = -%Ax = 8, g» = f lM ^' A - ) . (20) 

D A v c / d a c 2 d a V 6 y v 7 

Hence the cross section for two images with a splitting angle A9 > A9q is given by 

<7 = 7r^tf(A0- A0 O ) =ttCo^(M-M ) , (21) 

where $ is the step function, and M related with A9 can be solved from the equation (TSU]) 

The probability for a source at redshift z s undergoing a lensing event on account of the 
galaxies distribution from the source to the observer can be obtained by dividing the total 
lensing cross-section by the area A(z) of the lens plane, (Schneider et al. 1992) 

PZ S POO JJJ 

P= / — -(l+zfn(M,z)a(M,z) dMdz, (22) 
Jo Jo dz 

where is the proper distance from the observer to the lens. 

Inserting equation ([8]) and (l2Tj) into equation (122|) . we have for the SIS case (Li & 
Ostriker 2002) 

dP{> A9 ) _ dD L (D A D A LS 

107T p clit fi\l M {i + Z 



x J™ (^) 4 m - M )dM. (23) 

Differentiate this expression with respect to A9q, we can obtain the probability density for 
a source to have a double image with splitting angle A9 = A9o: 

d 2 P(> A9 ) 3 ^ dD L (D A D A S \ 2 f(M ,z) ( a v (M Q )\ 4 dM { 

1D7T p C rit,0"MU + Z 



dA9 dz — • r^,u--,« v - . ~, (iz \ D A J M \ c J dA9 ' 

dM n 3c 2 /—a \ 2 / 3 ni 



^M^ 3 . (24) 



dA# 16ttG V 6 

The proper distance D L and the angular-diameter distance D A from the redshift Z\ to 
22 are calculated via 

D L (ZU Z2) = I" TT-^TT , D A ( Zl , Z 2 ) = T^—T4h- ^ 

J Z1 (1 + z)H(z) 1 + z 2 J Zl H(z) 
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3.2. NFW Profile as a Lens 

The NFW profile, based on the N-body numerical simulation of cold dark matter, is a 
very important approach for understanding the formation of galaxies and clusters of galaxies 
(e.g., Zhao 1996; Hanyu & Habe 2001). Its mass density is given by (Navarro, Frenk & 
White, 1995, 1996, 1997) 

3 

PnfwM = , P ^ s , 2 , (26) 
ryr + r s ) 2 

where p s and r s are constants: 



PM^vir c\ 1 / 3M 
Ps = o T7~^ ' r s = — ( - A x — ' ( 27 ) 



I 1 / 3AY 

/(ci) ' s ci \A7rp M \ 

with /(ci) = ln(l + ci) — Ci/(1 + Ci). For the concentration parameter ci, we adopt the 
fitting formulae given by Bullock et al. (2001): c x = 9(1 + ^(M/l.S x lO 13 ^-^©)- 013 . 

Similar to the case of SIS, we define the position vector in the lens plane and the source 
plane as £ = xr s and t] = y r s Dg / , respectively. The surface mass density for the NFW 
profile is given by(Li & Ostriker 2002) 

E(x) = 2p s r s / (x 2 + z 2 )- 1 ' 2 ^ 2 + z 2 ) 1 ' 2 + l]- 2 dz (28) 
Jo 

Then the reduced lensing equation is 

y = x- p s ^- , (29) 

x 

where 



Ps = 4p s r s /S cr , 
1 D A 



s 



3C 



s(x)= I udu I (u 2 + z 2 ) 1/2 (u 2 + z 2 Y rz + l ~dz, (30) 



2\l/2 



-2 



The dimensionless parameter p s determines the size of the lensing cross-section a for a NFW 
halo to produce multiple images: larger p s , smaller a. The curve of y to x runs through 
the coordinate origin and has a extremum point central-symmetrically on each side, whose 
coordinates (x cr ,y cr ) are determined by dy/dx\ Xcr = and y cr = y{x cr ). Thus a single source 
with a certain y has multiple images when \y\ < y cr . Once more than two images are formed, 
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we shall only consider the splitting angle A9 between the two outside images. According to 
Li & Ostriker 2002, we shall neglect the variety of A9 caused only by the movement of y and 
get Ax (y) ~ Ax (y — 0) = 2x , where x is the positive root of y(x) = 0. Then the splitting 
angle AO is given by 

and the cross-section for forming multiple images with A9 > A# is 

a (> A9 , M, z) « ny 2 CI r 2 s {AO - A6 ) . (32) 

Fig. 3 shows the splitting angle AO as the function of M(lO 15 /i" 1 M ) in SIS and NFW 
cases for w = —0.5, w = —1.0, w = —1.5, and w(z) = —1.0 — z/(l + z). The source object is 
at z s = 1.5, and the lens object is at z = 0.3. In Fig. 4 we plot the A9 against the redshift of 
lens z for M = 0.01, 1.0, and 100, respectively. The source object is at z s = 1.5. The curves 
of w(z) and w = —1.0 almost overlap in both figures. We can see that the A6 produced by a 
NFW lens is more sensitive to the parameter w than that of an SIS lens, especially for small 
M. For SIS case, There are only quite small changes of A9 for our different selections of w 
in both Fig. 3 and Fig. 4. 

Fig. 5 gives the cross section a against the redshift of lens z for M(10 15 h~ 1 M ( 7 ) ) = 0.01, 
1.0, and 100, respectively. The source object is at z s = 1.5. For SIS case, There are not 
visible changes of A9 for our different selections of w. For NFW case, the cross section is 
more sensitive to the parameter w for smaller M and somewhat smaller z. 

Using equations (jSJ) , (122]) . and fl32l . we then get the differential lensing probability for 
the NFW case (Li & Ostriker 2002) 

dP(>A9 ) _ n . 3 dD L f°°f(M,z) 



dz 



7Tp C ru,o n M (l + zf ^ ^ l^p-ylr^{M - M )dM . (33) 
dz J M 



Differentiate this expression with respect to A6 , we can obtain the probability density for 
a source to have a double image with splitting angle A9 = A9 : 

d 2 P(> Ag ) 3 dDi/(JM 2 2 (JM 

dA9 dz = -P^ M (l + z) ^^^UMoK(M )— , 

dM Dj; ( , n/r ,dx l/3 + 0.13 N_1 



dA0 2 y r,(M a )— +Xo (M )r,(M a ) ^ 

dx g(x )/x dfx s /dM 

dM ~ 1- fi s (M )(g'(x )/x -g(x )/x 2 ) ' 
dfj, a 4r s p s / 0Q7+ 0.13c? 



dM M E cr V ' (l + Cl ) 2 /(ci) 



X 



(x 



1/2 



(34) 



4. DATA ANALYSIS AND NUMERICAL RESULTS 

The CLASS statistical sample provided a well-defined statistical sample with N = 8958 
sources. Totally iVj = 13 multiple image gravitational lenses have been discovered and all 
have image separations A8 < 3" (Browne et al. 2003). The data informations of the 13 
observed lens systems are not entire: the source redshift z s and lens redshift z are both 
unknown for 1 lens systems, only z s unknown for 4 lens systems and only z unknown for 1 
system. 



4.1. Basic Preparations 

The CLASS statistical sample use the flux density ratio q r of the multiple lensing images 
as a selection criterion of a sample: q r = < 10 (Chae et al.2002; Chen 2003a, 

2003b), where the //+ and //_ are the magnifications of two (outside) images, respectively. 
The magnification of a image is determined by \i = [f ^f] • The parameter q r reduces the 
lensing cross section a. For a NFW case, the influence of q r on a is very small and can be 
neglected. For an SIS lens, it needs to multiply the a by a factor (9/1 1) 2 . 

As the exact redshift distribution of the CLASS statistical sample is unknown, Chae et 
al.(2002) utilized a Gaussian model with mean redshift (z 8 ) = 1.27 given by Marlow et al. 
(2000) to describe the redshift distribution for the unlensed sources of the CLASS statistical 
sample. Denote this Gaussian model as g c (z s ), its distribution was explicitly plotted in 
Figure 5 of Chae (2003), which was obtained from describing the redshift distribution of the 
flat-spectrum sources as shown in that Figure 5. Note that as such an Gaussian model has 
a physical cut at the point z s = 0, it is no longer to have the standard form. Unlike the 
Gaussian model g c {z s ) by Chae (2003), here we shall take an alternative Gaussian model 
by directly fitting the redshift distribution of the subsample of CLASS statistical sample 
given by Marlow et al. (2000) instead of fitting the redshift distribution of the flat-spectrum 
sources in Chae (2003). Taking the general form of Gaussian model 



, , N. 

3(z ' } = 7!S exp 



[z s - a) 2 



2A 2 



(35) 



with N s being the normalization parameter J °° f{x)dx = 1, and requiring the mean value 
f °° xf(x)dx = 1.27 given by Marlow et al. (2000), we then only need to fit the remaining 
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one parameter. The best fit results are found to be 

N s = 1.6125; a = 0.4224; A = 1.3761 (36) 

For comparison, we also present some results based on two treating methods appearing 
in literatures by using the CLASS statistical sample: the redshift distributions of sources 
are the average redshift value d(z s ) = 5(z s — 1.27) (see, e.g., Li & Ostriker 2001) and 
f(z s ) = 0.204 + 0.2979,2, - 0.1121^ + 0.00 1584z s 3 (see, e.g., Sarbu et al. 2001). Fig. 6 gives 
curves of four models as a function of z s and the histogram of 27 CLASS subsample from 
Marlow et al. (2000). 

Before comparing with the CLASS statistical sample, we should consider the effect of 
magnification bias B, which causes the overrepresentation of the lensed objects in a flux- 
limited survey. The flux distribution of the CLASS statistical sample is well-described by 
N(f) oc (f/f ) v with rj = 2.07±0.02(1.97±0.14) for / > f (f < f ) and /„ = 30mJy (Chae 
et al.2002). The analysis process of magnification bias is determined by the lensing equation 
and following equations 



B 



JT«) 

A = ZAi 



(37) 



ydy dy/dx\i 

P(A) = ^L- ;p (y) = a y ] J p(y)dy = l, (38) 

where x and y are the position parameters on the lens plane, A is the total amplification 
of the multiply-imaged sources, A m is the minimum value of A and p(A) is the probability 
density for the amplification A, y max is the maximum value of a source' position on the 
lens plane, a is the normalization constant of probability function p(y): for the SIS case 
Umax = 1, « = 2 and for the NFW case , it can be calculated out numerically by using 
dy/dx\ Vmax = 0. Then for the SIS case, after considering the influence of flux density ratio 
q r , one needs only multiply the lensing probability by a constant factor 3.36. For the NFW 
case, B can be obtained from the numerical calculation. P Q b s (> A0) and dP ohs (> A9)/dA9 
are related to P by an integration 



and 



p(w) = P obs (> A0) = J j B dP{> d M) ^{z s )dzdz s , (39) 

dP ohs (> AO) f f d 2 P(> Ag) 
?H - dM =JJB — <p(z.)dzdz. , (40) 
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with <f(z s ) is the redshift distribution of sources. 

We shall compare the theoretical results of the SIS case and NFW case with the CLASS 
statistical sample. Fig. 7 shows the lensing probability P(> AO) as a function of the 
splitting angle AO for the source redshift distribution d(z s ), f(z s ), g c (z s ) and g(z s ). The 
thickest line in each panel is induced from the 13 observed lensing data. It is seen that the 
lensing probability P(> AO) of SIS case and NFW case are both sensitive to the parameter 
w, especially for NFW case. When parameter w increases, the values of P(> AO) for both 
SIS and NFW cases clearly increase in the whole concerned ragion of AO. Thus it is feasible 
to constrain the parameter w from the SGL splitting angle data. The figure shows that the 
SIS model can only reproduce the data curve at small AO < 1.5". When consider the rapid 
decline of P(> AO) from the data line at large AO, a combined mechanism of SIS and NFW 
model is needed to explain the whole experimental curve. Define a new model parameter M c 
as Li & Ostriker (2001): lenses with mass M < M c have the SIS profile, while lenses with 
mass M > M c have the NFW profile. Then the differential probability 

dP/dM = dP S is/dMd(M c — M) + dP NFW /dM§(M - M c ) 

where i? is the step function, d(x — y) — 1, if x > y and otherwise. Because the splitting 
angle AO is directly proportional to the mass M of lens halos, the contribution to large AO 
of SIS profile is depressed by M c . The lens data require a mass threshold M c ~ 1O 13 /i _1 M , 
which is consistent with the halo mass whose cooling time equals the age of the universe 
today. In this note we shall use such a two-model combined mechanism to calculate lensing 
probabilities. 

In Fig. 7, One can also find the influences of different source distributions on the 
lensing probability P(> AO): the patterns of the function curves are hardly changed, but 
the function values for the same AO slightly increase from d(z s ) to f(z s ). 



4.2. Constraint on w 

We now come to our main purpose i.e., utilizing the SGL splitting angle statistic data 
from strong gravitational lenses to constrain the equation-of-state parameter w of dark en- 
ergy. For that, we define the likelihood function as 

L(w) = (l-p(w)) N ^Hq t (w). (41) 

i=l 

p(w) represents the model-predicted lensing probabilities P(> 0.3") of a source with the 
redshift distribution <p(z s ) and can be calculated by using equations ff23l . fl33|) . and (139]) . 
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Here A9 > 0.3" is an observational selection criteria. qi{w) is the model-predicted differential 
lensing probabilities dP(> A9i)/dA9i of i th observed lens system with the splitting angle 
A9i and can be calculated by using equations ($3$, (13~4"|) . and ( |4"0"]) . To utilize the data 
informations adequately, we multiply the equation (J35|) by 5(z — z l ) or/and S ( z s — z\) (replace 
<p(z s )) f° r the i th observed lens system whose lens redshift z l or/and source redshift z\ is/are 
known. For the unknown z % or z l s , we just integrate it out. Compared with utilizing the curve 
of lensing probability P(> A9) as a function of the splitting angle A9 in the last subsection, 
which is introduced from lensed signals only, by using the likelihood function equation (jHJ), 
the unlensed signals are utilized and their influences are quite important due to the large 
exponential number (N — Nl) = 8945. 

Firstly we discuss the possible constraints on the model parameter M c and the constant 
w. Under the given cosmological parameters (f2 M ,/i,cr 8 ) = (0.24,0.73,0.74), Fig. 8 shows 
the 68% CL. and 95% CL. allowed regions from the CLASS statistical sample for the 
source redshift distribution g(z s ), g c (z s ), d(z s ) and f(z s ), respectively. The crosshairs in 
three panels mark the best-fit points (w, M c ) = (-0.89, 1.37), (-0.94, 1.36), (-1.4, 1.68) and 
(-0.73, 1.27) from left to right, where the unit of M c is lO^/i^M©. The source redshift 
distributions g(z s ) and g c (z s ) give nice constraints, while d(z s ), which is not proper redshift 
distributions, provides somewhat strange unexpected results. The redshift distribution f(z s ), 
which has a larger part of galaxies at high redshift, prefers a larger w and smaller M c . From 
now on we will only consider the g(z s ) case and g c {z s ) case, which are extracted from the 
subsample of the CLASS statistical sample (Marlow et al. 2000). Our best fit result of 
M c « 1.40 for both g(z s ) case and g c {z s ) case is larger than the value M c « 1.0 obtained by 
Li & Ostriker 2002. The 95% CL. allowed regions of parameter w for g(z s ) case and g c {z s ) 
case are from —0.18 to —1.45 and from —0.11 to —1.85, which are consistent with the ACDM 
cosmology. The Fig.2 in Chae (2007) shows a much negative result of the parameter w. Our 
consideration here is based on the specially selected cosmological parameters, and also we 
have used the two-model combined mechanism, namely the utilization of model parameter 
M c . As a consequence, our results avoid the large absolute value of parameter w. 

In Fig. 9, we show the constraint for the parameters (wq, w a ) appearing in a time- varying 
equation-of-state w(z) = w +w a z/(l + z) under the given cosmological parameters (fi^, h) = 
(0.24,0.73) for both g(z s ) and g c (z s ) cases, and with two different values of a 8 = 0.74 and 
erg = 0.90. The crosshairs mark the best-fit points (M c ;wo,w a ) = (1.36; —0.92, —1.31) for 
g(z s ) case and (M c ;wo,w a ) = (1.38; —0.89, —1.21) for g c (z s ) case when erg = 0.74, and 
{M c ;w ,w a ) = (1.56; -0.81, -2.5) for g{z s ) case and {M c ;w ,w a ) = (1.54; -0.83, -2.22) 
for g c (z s ) case when erg = 0.90. The parameter erg has significant influences on the mass 
power spectrum and the number density of dark halos, and our results show that the best 
fit (wq, w a ) are changed for the different selections of the parameter as'- when as increases 
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from cr 8 = 0.74 to ag = 0.9, the best fit parameters wo increase moderately and the best fit 
parameters w a have a sizable decrease for both the g(z s ) and the g c (z s ) cases. 

After marginalizing the cosmological parameters (Qm, h, 0g) and the critical mass pa- 
rameters M c by Monte Carlo method, we obtain the constraint on (u> , w a ) in Fig. 10. 
For the three cosmological parameters, we assume the Gaussian prior distributions induced 
from the results of WMAP Three Year Data: (Qm ± cr m , h ± ah, erg ± a CT8 ) = (0.238 ± 
0.019,0.73 ± 0.03,0.74 ± 0.06). For the model parameter M c , we integrate it from 1.0 to 
4.0. The crosshairs mark the best-fit point (w ,w a ) = (—0.88,-1.55) for g(z s ) case and 
(wo,w a ) = (—0.91, —1.60) for g c (z s ) case. At the 95% C.L., our fitting results are consistent 
with that of Barger et al. (2006), and the SGL splitting angle statistic with the source 
redshift distribution g(z s ) and g c (z s ) gives somewhat more negative values for the parameter 



5. CONCLUSIONS 

From the above analyzes, we have shown how the SGL splitting angle statistic can be 
used to quantitatively constrain the equation-of-state parameter w of dark energy. Though 
due to the limited space-time, the difference of the parameter w has few influences on the 
splitting angle A8 and lensing cross section er, while through the comoving number density of 
dark halos as sources and lenses described by Press-Schechter theory, dark energy can affect 
the efficiency with which dark-matter concentrations produce strong lensing signals. With 
a two model combined mechanism of dark halo density profile, which introduces a model 
parameter M c , we have carefully investigated the constraints on constant w and time-varying 
w(z) = w +w a z/(l+z). We find the best fit value M c ?a 1.4 for both w and w(z), such a value 
is larger than the value M c « 1 obtained by Li & Ostriker 2002. This is mainly because in our 
analyzes both the lensed and unlensed signals have been utilized in the likelihood function 
equation f|4"Tl) . while in the analyzes by Li & Ostriker 2002, only the lensed signals were be 
used. The transition from SIS to NFW characterized by the parameter M c is also motivated 
by the process of baryonic cooling (e.g., Kochanek & White 2001), where the parameter M c 
was introduced to divide the cooled(SIS) and uncooled(NFW) halos. The estimated value 
by Kochanek & White (2001) for h = 0.67 is M c w 1 x 10 13 M stm for the model without a 
bulge, which is also smaller than our fitting result. There are several differences between 
our calculations and theirs. Firstly, we use a larger Hubble constant h = 0.73 > h = 0.67 
and in our Figure 9, we show that for a larger h, we have a larger fit M c . Secondly, the SIS 
model used by us has a larger relative lensing cross section than the exponential disk used by 
Kochanek & White (2001), which gives a larger M c . Thirdly, the ratio of the cluster to the 



-17- 



background density A vir = ptaio/ Pm > 150 used by us is much larger than the value A ps 100 
used by Kochanek & White (2001), which could also have some influences on the results. 
Nevertheless, all the fitting results for M c are consistent with appropriate considerations. 

With the given cosmological parameters (£lM,h,a 8 ) — (0.24,0.73,0.74), we have com- 
pared the results of constant w corresponding to four kinds of source redshift distribu- 
tions. It has been shown that d(z s ) is not suitable for the SGL data analysis. For the 
redshift distributions of normalized Gaussian-type model g(z s ) and Gaussian model g°(z s ), 
the fitting results are (w, M c ) = (-0.89^, 1.37±g;||) and (w, M c ) = (-0.94+°;^, 1.36+°;^) 
respectively, and the fitting results for the constant w are consistent with the ACDM at 
95% C.L.. For the time-varying w(z), we have firstly investigated the influence of as with 
the redshift distributions g{z s ) and g c (z s ) and found that the fitting results of the double 
parameters (w ,w a ) are changed when a 8 increases from a 8 = 0.74 to a 8 = 0.9. With 
the above given cosmological parameters, the best fitting results for the g(z s ) case are 
(M c ;w ,w a ) = (1.36; -0.92, -1.31) for a 8 = 0.74 and (M c ;w ,w a ) = (1.56; -0.81, -2.5) 
for og = 0.9; and for g c (z s ) case, the best fit results are (M c ; wo,w a ) = (1.38; —0.89, —1.21) 
for a 8 = 0.74 and (M c ; w ,w a ) = (1.54; -0.83, -2.22) for a 8 = 0.9. 

After marginalizing our likelihood functions over the cosmological parameters (Qm, h, <r 8 ) 
(by using the prior probabilities induced from the WMAP Three Year data) and the model 
parameter M c , we have obtained a reliable constraint on the parameters (wo,w a ). Within 
the allowed uncertainties, the results for wq are consistent with the constraints obtained 
from the Type la supernovae data (Barger, Guarnaccia & Marfatia 2006). Our fit results 
are (w ,w a ) = (-0.88+?;^, -1.551^) for g(z s ) case and (w ,w a ) = (-0.911?;^, -1.60±|;|?) 
for g c (z s ) case. It is noticed that the best fitting results based on the SGL splitting angle 
statistic favor negative values for the parameter w a , which differs from the best fitting values 
obtained based on the Type la supernovae data, where the best fitting results favor positive 
values for the parameter w a (Barger, Guarnaccia & Marfatia 2006). A combining constraint 
is interesting and will be investigated elsewhere. 

In conclusion, the quantitative investigation has shown that the SGL splitting angle 
statistic can lead to a consistent constraint on the constant w and the double parameters 
(w , w a ) of the time- varying dark energy equation of state w(z) = w +w a z/(l+z). Especially 
for the allowed range of parameters (wo,w a ), the SGL splitting angle statistic does give an 
interesting bound. It can be seen from Fig.8 to Fig. 10 that the normalized Gaussian-type 
source redshift distribution g(z s ) leads to the most stringent constraints. Though it does 
not yet allow to obtain a more accurate constraint, while it can provide a complementarity 
to other constraints from Supernovae, cosmic microwave background, weak lensing. 
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Fig. 8. — 68% C.L. and 95% C.L. allowed regions from the CLASS statistical sample for the source 
redshift distribution g(z s ), g c (z s ), d{z s ) and f(z s ), respectively. The crosshairs in three panels mark 
the best-fit points (w,M c ) = (-0.89,1-37), (-0.94,1.36), (-1.4,1.68) and (-0.73,1.27) from left 
to right. 
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Fig. 9. — 68% C.L. and 95% C.L. allowed regions from the CLASS statistical sample for the 
f{z s ) redshift distribution with different ag = 0.74 and 0.9, respectively. The crosshairs mark 
the best-fit points: for the g(z s ) case are (M c ;wo,w a ) = (1.36; —0.92, —1.31) for ag = 0.74 and 
(M c ;wo,w a ) = (1.56; —0.81, —2.5) for ag = 0.9; and for g c (z s ) case, the best fit results are 
(M c ;w ,w a ) = (1.38; -0.89, -1.21) for cr 8 = 0.74 and (M c ;w ,w a ) = (1.54; -0.83, -2.22) for 
a 8 = 0.9; 
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Fig. 10. — 68% C.L. and 95% C.L. allowed regions from the CLASS statistical sample for the source 
redshift distribution f(z s ). The crosshirs mark the best-fit points (wo,w a ) = (—0.88,-1.55) for 
g(z s ) case and (wo,w a ) = (—0.91,-1.60) for g c (z s ) case . The cosmological parameters (£Im, h,as) 
and the model parameter M c have been marginalized 



